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We study the two and four dimensional Nishimori mul- 
ticritical point via high temperature expansions for the ±J 
distribution, random-bond, Ising model. In 2d we estimate 
the the critical exponents along the Nishimori line to be 
7 = 2.37 ± 0.05, v = 1.32 ± 0.08. These, and earlier 3d esti- 
mates 7 = 1.80 ± 0.15, v — 0.85 ± 0.08 are remarkably close 
to the critical exponents for percolation, which are known to 
be 7 = 43/18, v = 4/3 in d = 2 and 7 = 1.805 ± 0.02 and 
v = 0.875 ± 0.008 in d = 3. However, the estimated 4d Nishi- 
mori exponents 7 = 1.80 ± 0.15, v = 1.0 ± 0.1, are quite 
distinct from the Ad percolation results 7 = 1.435 ± 0.015, 
v = 0.678 ± 0.05. 



In recent years there has been much interest in the 
study of critical phenomena in quenched-random, two- 
dimensional, thermodynamic systems. However, with the 
exception of percolation, for which various critical pa- 
rameters are known exactly, other random hxed points 
are not fully understood. Such random critical phenom- 
ena are of interest from a theoretical point of view, and 
also from an experimental point of view. Among notable 
experimental systems showing two-dimensional random 
critical phenomena are pleateau transitions in quantum 
hall systems |jj and Bose-glass transitions in dirty super- 
fluids and superconductors 

Perhaps the simplest theoretical model with quenched 
randomness is the random-bond Ising model. In d = 2 
a lot is known about weak randomness ||. The case 
where randomness has the most dramatic influence on 
thermodynamic properties is that of a symmetric dis- 
tribution of bonds, that is one where there are roughly 
equal tendencies for ferromagnetic and antiferromagnetic 
ordering. In this case the system may only have a long- 
ranged spin-glass phase at low temperatures. There is 
considerable numerical evidence that in d = 2 there is no 
finite temperature spin-glass phase Q. The Nishimori 
manifold separates the region in parameter space where 
ferromagnetic ( or antiferromagnetic ) correlations are 
stronger from those where spin-glass correlations domi- 
nate. In certain random-bond Ising models many exact 
results can be obtained along this special manifold ||. 
A Nishimori multicritical point can exist even in the ab- 
sence of a finite temperature spin-glass transition and has 



been studied by renormalization group Q and various nu- 
merical methods . By now the existence of the critical 
point and its location are reasonably well established , 
although to our knowledge no reliable estimates of the 
critical exponents exist. 

Here we study this model by high temperature expan- 
sions, estimating the critical point and the various crit- 
ical exponents. Our estimates for the critical tempera- 
tures are consistent with previous ones. Our interesting 
result is that the critical exponents 7, and v are remark- 
ably close to that of percolation. This turns out to be 
the case in three dimensions also. To see if such a trend 
continues with dimensionality we study the four dimen- 
sional Nishimori multicritical point. There the critical 
exponents are clearly distinct from percolation. Towards 
the end of this paper we speculate on the relevance of the 
percolation fixed point to the present problem. 

We consider the Hamiltonian 

<i,j> 

with the Jij independent quenched random variable with 
distribution 

P{J l0 ) = pSiJij - J) + (1 - p)5(Jij + J) (2) 

The Nishimori line in the parameter space of temperature 
T and ferromagnetic bond concentration p is given by 

v = 2p - 1 (3) 

with v = tanh J/kT. It is most convenient to directly 
develop the expansions along the Nishimori line In 
this case the expansion variable is w = v 2 . We define 
susceptibilities Xm.n by the relations 

x fflltl =iE[< s < s i>T ( 4 ) 

Here angular brackets represent thermal averaging and 
the square brackets an averaging with respect to the bond 
distribution. We carry out expansions to 19th order in 
2d and 15th order in Ad for Xi,\ an d X2,2 using the star- 
graph method |9). We note that along the Nishimori 
line the ferromagnetic susceptibility, xi,i exactly equals 
the spin-glass susceptibility X2,i- In order to study the 
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crossover exponents we also calculate the series for Xv = 
i'c'X2,i/<9t>.The expansion coefficients are given in table 
1. 

We shall assume that the quantities x = X2,i, x' = 
X2,2 and Xv, become singular at a critical point w c with 
exponents 7, 7' and 7" respectively. Assuming standard 
scaling, the exponents for the divergence of the different 
series can be related to the critical exponents v, rj and (f> 
as §: 

7 = (2 - rfjv, 7' = (d - 4 - 277)2/, 7" = 7 + <j) 

Here d is the dimensionality of the system. 

We analyze the series based on the expectation that 
near the critical point the susceptibilities have the form 

X oc (w c — w) r (l + a(w c — w) Al + b(w — w c ) + •••)• 

We estimate the location of the critical point, w c , the 
dominant exponent 7, and the correction-to-scaling ex- 
ponent Ai. The value of Ai may be biassed by the pres- 
ence of even higher order correction terms, however fit- 
ting to this form should ensure reliable evaluation of the 
critical point and the dominant exponent. Our analysis 
is carried out with no prior assumptions regarding expo- 
nent values, and was done by one of us without knowing 
prior literature values for the critical parameters or scal- 
ing relations between the exponents. We have studied 
the series with two methods, commonly known as Ml 
and M2 |L0| . They are based on suitable transformations 
of the series, and Pade approximants for the transformed 
series. 

Ml: In this method of analysis we study the logarithmic 
derivative of 

B(w) = hH(w) - (w c - w)^^-, 

aw 

The dominant singularity is a pole at w — w c with a 
residue (h — 1). 

We implement method Ml as follows: for a given value 
of w c we obtain Ai versus input h for many central 
and high Pade approximants, and we choose the triplet 
w c , h, Ai where all Pades yield as nearly as possible iden- 
tical values of h. 

M2: In the second method we first transform the series 
in w into a series in the variable y, where 

y = l-(l-w/w c ) A \ 

and then take Pade approximants to 

G(v)=A 1 ( y -l)±Jn(H{w)), 
dy 

which should converge to —h. Here we plot graphs of h 
versus the input Ai for different values of w c and again 
choose the triplet w c , h, Ai, where all Pades converge to 
the same point. 



We found good convergence for all three d = 2 series, 
with M2 graphs being better converged everywhere, and 
Ml giving consistent results. The x an d x' series behaved 
better than the Xv series. For brevity, we only show some 
representative plots from which we have deduced our es- 
timates of critical parameters. In Figure 1 we present two 
three dimensional graphs from the M2 analysis for the x 
series. In Figure 1(a) we show the three dimensional ver- 
sion on a fairly coarse temperature scale, in Figure 1(b) 
we show a finer scale. On the coarse scale we can see that 
at the trial w c values of 0.525 and 0.550, convergence to a 
region of clear intersections is much poorer than at 0.575. 
Similarly, although the different approximants come to- 
gether at the background plane with w c = 0.625, this 
type of almost flat graph with the asymptotic conver- 
gence at very high Ai values is indicative of behaviour 
that does not give correct critical behaviour in test sys- 
tems or exactly solved models. (It is often seen near trial 
critical points that give exponent values that violate hy- 
perscaling.) The fine scale shown in the enlargement in 
Figure 1(b) shows us a set of graphs which mostly sat- 
isfy the considerations required of an intersection region, 
with the best of all being the central plane. Thus from 
the M2 analysis, the best w c estimate is 0.596 ± 0.008. 
This implies p c = 0.886 ± 0.003, T c /J = 0.975 ± 0.006, 
which are consistent with previous estimates . In Fig- 
ure 2(a) the central slice at w c — 0.596 is shown. From 
this we conclude an exponent estimate of 7 = 2.37±0.05. 
The Ml analysis is consistent with these values. In the 
M2 analysis for the x' series, convergence was again op- 
timal at w c — 0.596. The exponent, 7', is deduced to 
be 2.11 ± 0.07. Via scaling this gives v = 1.32 ± 0.08. 
We analyze the Xv series in two ways, first by consider- 
ing Xv/w arL< i second by studying the series for dxv/dw. 
From these analyses we conclude 7 + = 3.0±0.3. 

In four dimensions, the M2 analysis of the x series gives 
best convergence at w c = 0.1764, with 7 = 1.9. From Ml 
a slightly lower w c = 0.176 value and a correspondingly 
lower 7 = 1.8 seems optimal. For the %' series optimal 
convergence is at w c — 0.176 with the central values of 
the exponent 7' of 0.41 from Ml and 0.40 from M2. Since 
three of the four analysis support the lower value of w c — 
0.176 we use this for our final estimates. We quote overall 
values of w c = 0.176 ± 0.001 and 7 = 1.80 ± 0.15, 7' = 
0.40 ± 0.03 By scaling this gives us v = 1.0 ± 0.1. 

As stated in the introduction, these estimated critical 
parameters are remarkably close to percolation [Tlfl in 
d = 2. The numbers for percolation are 7 = 43/18 and 
v = 4/3. Furthermore, in d = 3 the Nishimori exponents 
were found to be 7 = 1.80 ±0.15 and v = 0.85 ±0.08 [§. 
These numbers are also confirmed by our present analy- 
sis. The d — 3 percolation exponents are 7 = 1.805 ±0.02 
and v = 0.875 ±0.008 0. Thus, in d = 3 also the Nishi- 
mori exponents are very close to percolation. However, 
in d = 4 the percolation exponents are 7 = 1.435 ± 0.015 
and v — 0.678±0.05, which are clearly distinct from those 
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found here. We note also that the exponent v along the 
Nishimori line appears to be non-monotonic as a function 
of dimensionality. 

One can conclude that in d = 2 and 3 the Nishimori 
critical behavior is consistent with the universality class 
for percolation. However, this is not so in d = 4. The 
possibility that the closeness of the Nishimori exponents 
to percolation in d = 2 and 3 is purely accidental can- 
not be ruled out. However, the following considerations 
suggest a possible connection. If we consider a bond dis- 
tribution of the form: 

P{Jij) = pS(Jij - J) + (1 - p)S(Jij + J) + cS^) 

That is, the bonds are allowed to take values ± J as well 
as zero, then in the generalized parameter space of p, 
c and temperature T, the Nishimori manifold is a two- 
dimensional plane. For p — ( or unity), this plane 
reduces to the T = dilution axis and thus contains 
the percolation fixed point |12|j . However, it is generally 
believed that T is always an unstable direction for perco- 
lation and, hence, finite temperature Nishimori criticality 
should have different exponents. 

Secondly, Nishimori has argued (Tj| that the spin-glass 
ferromagnetic transition is a geometry-induced phase 
transition ( as opposed to a thermal transition ) , which is 
also true for percolation. However, if this led to the iden- 
tification of the Nishimori fixed point with percolation it 
should be true independent of dimensions. However, our 
results in d — 4 contradict this. Furthermore, the epsil- 
lon expansions (around d = 6) for the exponents at the 
Nishimori multicritical point are different from percola- 
tion §. 

Finally, a very different way in which this model is of 
significant interest, is through the mapping between the 
2d Ising model and free fermions in 1 + 1-dimension, and 
the connection between the latter and the pleateau tran- 
sitions in the quantum hall effect. The Chalker model for 
the pleateau transitions in the quantum hall effect [ ^4|Jl5| 
can thus be mapped onto random-bond Ising models. 
However, the ± J model studied here does not have the 
correct symmetries for the quantum-hall problem [JL6| . 
It is well known that percolation occurs in one limit of 
the quantum-hall systems, when the disorder potential is 
slowly varying in space |l7| . However, numerical studies 
of the Chalker models lead to exponents clearly different 
from percolation |l4| . 

Thus, there is no compelling theoretical reason why 
the Nishimori multicritical point should be in the uni- 
versality class of percolation. One possible explanation 
for our findings could be that in low-dimensions, where 
the multicritical point occurs at very low temperatures, 
there are crossover effects which produce effective expo- 
nents close to percolation. These issues deserve further 
attention. 
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FIG. 1. Plots from the M2 analysis of the d — 2 \ series 
(a) on a coarse scale and (b) on a fine scale. 



FIG. 2. The central slice for the M2 analysis in Figure 1 
with w c — 0.596. 



TABLE I. Expansion coefficients for the susceptibilities in 
d = 2 ( See Eq. 4) 
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_n X20 X2£ Xv_ 

Oil 

14 8 

2 12 4 48 

3 36 216 

4 76 36 512 

5 196 -32 1640 

6 316 236 1856 

7 884 -464 9208 

8 780 1988 -5824 

9 3684 -5072 57576 

10 396 17076 -109264 

11 22740 -50432 680376 

12 -22596 164108 -1547376 

13 188420 -500496 8163624 

14 -331108 1604572 -21618432 

15 1517396 -5042160 84085560 

16 -4509268 16221028 -311253632 

17 15654148 -52336864 1071437960 

18 -56714548 170687620 -4263416944 

19 183041524 -561493296 14787979576 



TABLE II. Expansion coefficients for the susceptibilities 
in d = 4 ( See Eq. 4) 



J}_ Xxi X2£ Xv_ 

Oil 

1 8 16 

2 56 8 224 

3 392 2352 

4 2552 200 19840 

5 16904 -192 162512 

6 105944 6584 1179840 

7 679784 -12384 8736880 

8 4158200 234824 58846592 

9 26120392 -649056 412368720 

10 157020984 8748712 2651405536 

11 974362408 -32109952 18054488432 

12 5783009304 342786296 112459651552 

13 35661616648 -1523180000 755621878608 

14 209506120728 14008147224 4590427798720 

15 1289118273320 -70814307872 30721400183024 
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